Quantum-Critical transport at a semimetal-to-insulator transition on the honeycomb 

lattice 



Lars Fritz 

Institut fiir Theoretische Physik, Universitdt Koln, Ziilpicher Strafle 77, 50937 Koln, Germany 

(Dated: January 26, 2012) 

In this paper we study transport properties of electrons on the two-dimensional honeycomb lattice. 
We consider a half-filled system in the vicinity of a symmetry-breaking transition from a semimetal- 
lic phase towards an insulating phase with either charge density or spin density wave order. The 
effect of either order is to break the sublattice inversion symmetry which induces a finite gap for 
the electronic single-particle excitations. Phenomenologically, such a scenario is described in the 
framework of a Gross-Neveu theory. We analyze two related formulations of the model by means 
of (i) a controlled renormalization group calculation and (ii) the large-N method, both of which in 
combination with a Boltzmann transport equation. We determine the quantum-critical conductiv- 
ity and also discuss crossover behavior from quantum critical behavior into the insulating and/or 
the semimetallic phases. We find that at asymptotically low temperatures the quantum-critical 
conductivity is given by a temperature independent universal number. Over a large temperature 
window the temperature independent quantum critical conductivity is masked by a logarithmically 
temperature dependent contribution due to the marginally irrelevant long-range Coulomb interac- 
tion. We discuss possible origins of this peculiarity in the two complementary formulations of the 
model. Furthermore, we consider possible relations of our findings to recent experiments, with a 
special emphasis on the quantum-critical-to-insulator crossover. We find that our results are in re- 
markably good qualitative and quantitative agreement with a recent analysis of the data sets under 
the hypothesis of an underlying gap in the single-particle spectrum. 



I. INTRODUCTION 

Two dimensional electronic systems on a honeycomb 
lattice have seen tremendous activity over the last couple 
of years. This development was triggered by the exper- 
imental isolation of single layer graphene which comes 
with a variety of fascinating and unusual propertieJ^. 

Commonly, graphene is considered an excellent con- 
ductor with extremely high mobilities. The effective 
low-energy theory describing the system is given in terms 
of massless Dirac fermions moving at the Fermi velocity 
vf instead of the speed of light c « 3^). Most 
of the experimental findings in graphene to date can be 
explained in terms of weakly interacting massless Dirac 
quasiparticleS. A most prominent exception is pro- 
vided by the recently observed fractional quantum Hall 
effecPa. 

Despite this predominant lack of signatures of correla- 
tion effects, many authors have analyzed the role of local 
and long-range Coulomb interaction and possible inter- 
action driven insulating phaseS^i^. 

In this paper we study transport in the vicinity of 
a semimetal-to-insulator quantum phase transition, in 
which the opening of the single-particle gap in the insu- 
lating phase is driven by spontaneous sublattice symme- 
try breaking. The sublattice inversion symmetry break- 
ing is driven by collective instabilities of the charge den- 
sity or spin density wave type. Note that this scenario 
does not describe the Mott transition in the traditional 
sense due to the additional symmetry breakinj^. In that 
sense the treatment presented here is closer to the exci- 
tonic insulator obtained due to long-range Coulomb in- 
teraction and chiral symmetry breakinglSIlH than it is to 



a pure Hubbard model. Within this paper we analyze 
the simplest possible transition from a semimetal to an 
insulator described by a Hubbard like model and discard 
the possible existence of interesting intermediate phases 
such as a spin liquid phas^^^. 

Transport properties at quantum-criticality have tradi- 
tionally been a very active branch of research in strongly 
correlated matter. In two dimensions microscopic calcu- 
lations have been performed with the quantum-critical 
conductivity being given by a temperature independent 
universal numbep22H23l_ This has to be expected for di- 
mensional reasons if the critical point of the field theory 
has relativistic invariance and is located at finite interac- 
tion strength. 

In the context of clean intrinsic graphene, which shares 
a number of characteristics with a quantum-critical sys- 
tenJ^, the critical conductivity was shown to diverge in 
logarithmic manner upon lowering temperature^'^. This 
behavior is rooted in the renormalization group fixed 
point of the system lying at zero interaction, with the 
interaction parameter describing long-range Coulomb in- 
teraction being marginally irrelevant, thus flowing to zero 
logarithmically upon lowering the energy scale. 

In our problem, the insulating phase is characterized 
by the presence of collective order of charge density 
(CDW) or spin density wave (SDW) type^o. The fi- 
nite order parameter acts like a mass term for the Dirac 
fermions and consequently opens a gap in the electronic 
excitation spectrum. Both sorts of long-range order 
break sublattice inversion (Z2) symmetry (convention- 
ally referred to as chiral symmetry), which is the driving 
force behind the finite electronic gap. SDW additionally 
breaks spin rotation symmetry (SU(2)). 
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Within this paper we choose a simple field theoretic 
formulation of the Gross-Neveu typ^^Siili^ which captures 
the collective instability and the opening of the electronic 
quasiparticle gap. We consider two related formulations: 
(i) a Landau- Ginzburg type bosonic order parameter the- 
ory coupled to Dirac fermions via a Yukawa-type coupling 
(henceforth referred to as model I) and (ii) a locally in- 
teracting theory of Dirac fermions (henceforth referred 
to as model II). In both cases there is a quantum phase 
transition of the semimetal-to-insulator type, which is 
second order and described by an interacting fixed point. 
In both cases an appropriate order parameter can eas- 
ily be defined, which is zero in the semimetallic phase 
and finite in the insulating phase. Thinking in terms 
of the Landau-Ginzburg description of phase transitions, 
the order parameter symmetry is 0{N) with = 1 (Ising 
or Z2) for CDW and TV = 3 (Heisenberg) for SDW. 

Microscopically, local interactions can trigger these 
instabilities: on-site Coulomb repulsion favors SDW, 
whereas repulsion between adjacent sites favors CDW. 
Fig. [1] sketches a phase diagram with the vertical axis 
being temperature and the horizontal axis the local in- 
teraction parameter. The strength of local repulsion, 
denoted U, serves as tuning parameter for the quan- 
tum phase transition. In the case of a quantum phase 
transition towards CDW, U denotes repulsive interac- 
tion between adjacent sites, whereas for the case of a 
transition towards SDW the repulsion is on-site. Our 
analysis proceeds along the lines of a combination of a 
semiclassical Boltzmann equation combined with renor- 
malization group a rguments (m odel I) and a large-N ex- 
pansion (model II)22H23l25|28lE0] Throughout the paper 
all mathematical expressions are explicitly shown for the 
CDW case. The SDW expressions differ in combinatorial 
factors due to the difference in order parameter symme- 
try, but not in their structure. However, qualitative dif- 
ferences occur at finite temperature due to the absence of 
long-range order in the SDW case. In the CDW case a fi- 
nite gap in the electronic excitation spectrum is stable at 
finite temperatures, whereas there is no hard gap for the 
SDW. This is discussed in great detail in the conclusions. 



A. Overview of the results 

We calculated the quantum-critical minimal conduc- 
tivity of graphene at the semimetal-to-insulator transi- 
tion using two different phenomenological models. Fur- 
thermore, we estimated crossover functions for the con- 
ductivity which obtain upon lowering temperature and 
entering the semimetal or respectively the insulating 
phase. Most interestingly, we estimate the crossover 
function for the quantum-critical-to-insulator crossover 
and compare it to experimental data in graphene. 

Considering two slightly distinct models we find seem- 
ingly contrasting transport properties: 

(1) Within model I in contrast to other relativistically 
invariant critical points at finite interaction strength in 



two dimensional systems, the d.c. conductivity is not 
independent of temperature, but instead seemingly di- 
verges upon lowering temperaturd^fiHSS] xhis surprising 
result obtains because of a conspiracy of matrix elements 
and kinematical constraints which prohibits electronic 
current relaxation from the Yukawa coupling to all or- 
ders in perturbation theory. This statement, however, is 
only strictly true if the electrons and bosons are treated 
as sharp quasiparticles at all times. It turns out that the 
minimal conductivity (at charge neutrality) is eventually 
determined by a marginally irrelevant operator, namely 
long-range Coulomb interaction. The critical transport 
to leading order in temperature turns out to be identical 
to that of a gas of hot Dirac electr ons inte racting solely 
via long-range Coulomb interactiorP^^HllHlI], 

(2) Within model II we find a universal temperature 
independent universal conductivity as has to be expected 
from dimensional reasoning. However, we find that the 
prefactor of the inverse relaxation time is extremely small 
compared to the one associated with the marginally ir- 
relevant long-range Coulomb interaction. In contrast to 
model I, the scattering from the CDW (SDW) order pa- 
rameter fiuctuations is not completely inactive, but still 
extremely small. A simple order of magnitude estimate 
suggests that for all experimental purposes the universal 
quantum-critical conductivity is masked by the current 
relaxation due to long-range Coulomb interaction. 

We comment on the seeming discrepancy between the 
transport properties of model I and model II and argue 
that the two different pictures are actually compatible 
and describe the same basic physics, with model I over- 
estimating the kinematical blocking. 

On a more technical note another interesting result ob- 
tains within model I, in which bosons and fermions are 
directly coupled, namely the absence of boson-drag ef- 
fects for the coupled Boltzmann equations. Boson-drag 
constitutes a serious complication in obtaining trans- 
port coefficients in field theories of electrons coupled to 
bosons, where the bosons are an effective degree of free- 
dom, whose dynamics itself is governed by the underly- 
ing electrons. An often employed approximation for the 
calculation of electronic transport properties in such sys- 
tems is that the bosonic system is assumed to be equili- 
brated on timescales relevant for electronic transport and 
drag effects are neglected. This approximation turns out 
to be exact in our case. 

Finally, the crossover behavior of the conductivity 
for the the quantum-critical-to-insulator crossover is dis- 
cussed for both forms of collective ordering, namely the 
CDW and SDW. As we argued before, in the case of 
CDW an Ising degree of freedom condenses, which is al- 
lowed at finite temperature. Thus the electronic degrees 
of freedom have a hard gap, which entails an exponen- 
tially suppressed conductivity. In the case of the SDW, 
on the other hand, a real quasiparticle gap only opens 
at temperature T — and finite temperature behavior is 
not governed by a hard gap, but rather by a "pseudogap" , 
leading to a power-law suppression of the conductivity. 
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FIG. 1: Schematic phase diagram. The dashed lines denote 
crossover lines, whereas full lines denote phase transitions. 
The different regions and the crossovers are discussed in the 
main text. The insulating region III differs for the scenarios 
CDW and SDW. Whereas for CDW there is a finite tem- 
perature phase transition of the Ising type, there is no finite 
temperature transition for the SDW due to the Hohenberg- 
Mermin- Wagner theorem with forbids the breaking of a con- 
tinuous symmetry at finite temperature in two dimensions. 



of coupled transport equations neglecting the effect of 
" boson-drag" is exact to linear order in the applied elec- 
tric field and proceed to analyze the Boltzmann equation 
for the massless Dirac particles. We find that inelastic 
scattering due to the Yukawa-coupling is kinematically 
forbidden (see Sec. IIIA2I to a// orders in perturbation 
theory. This results in the absence of a temperature in- 
dependent minimal d.c. conductivity. In Sec. |III B| we 
solve the Boltzmann equation of model II. In contrast to 
model I we find that there is a finite universal conductiv- 
ity which does not depend upon temperature. However, 
we find that the prefactor of the inverse scattering time 



is tiny. This leads us to Sec. Ill C where we compare the 
results of in model I with those in model II and argue in 
which sense the two results are compatible. In Sec. |IV| we 
discuss the minimal conductivity in the quantum-critical 
regime and make an estimate of crossover temperatures. 
In a final section (Sec|v]) we discuss crossover behaviors 
which obtain upon lowering of temperature. Further- 
more, we discuss the fundamental difference between the 
CDW and SDW cases with a special eye on pseudogap 
behavior in the SDW case. A discussion of experimental 
conditions and reference to puzzling results of recent ex- 
periments on suspended high mobility graphene samples 
concludes Sec. Ivl 



Both scenarios allow to determine crossover curves 
which show remarkable similarities to experiments per- 
formed on suspended high-mobility graphene sample^^. 
It appears that they are qualitatively and quantitatively 
widely compatible with a recent analysis of the data sets 
under the hypothesis of an underlying electronic quasi- 
particle gap and might thus be a natural starting point 
for the identification of the fitting parameters discussed 
in Ref. [331 



B. Outline 

In SeclUwe introduce the scenario for the semimetal- 
to-insulator transition and motivate two different phe- 
nomcnological models, called model I and model II. Dur- 
ing the course of the paper we will switch between model 
I and model II, since different aspects of the problem are 
more conveniently discussed in one or the other formula- 
tion. 

We review the renormalization group treatment of 
model I following Ref. [10] and also the large-N treatment 
of model II following Ref. We furthermore stress some 
of the main physical properties of the models. In Sec. |III| 
we introduce the general formalism of the Boltzmann ap- 
proach. We apply the formalism to both models. In 
Sec. [Ill A] we consider model I and solve a system of cou- 
pled Boltzmann equations for the fermionic and bosonic 
degrees of freedom. The bosonic degrees of freedom cap- 
ture the physics of the collective instability and are even- 
tually responsible for the opening of the electronic quasi- 
particle gap. We explicitly show that in this system 



II. SPONTANEOUS INVERSION SYMMETRY 
BREAKING IN THE FRAMEWORK OF TWO 
DIFFERENT MODELS 

It is well known that a tight binding description of 
non-interacting electrons on the honeycomb lattice in its 
low-energy version can be cast in the form of the Dirac 
theory for massless fermions. The Hamiltonian for the 
electronic part readJ^ 

H^ivF J d^a;^^,«:7p9^*<T,«: , (1) 

with vf being the Fermi velocity. We explicitly consider 
half- filling, i.e., chemical potential fi = throughout the 
paper. Double indices are henceforth summed over unless 
stated otherwise. In the above notation we take ^'o-.k to 
be a spinor whose two components are associated with 
the two sublattices A and B. The index a denotes the 
spin degree of freedom whereas k denotes the valley {K 
and AT' point in the Brillouin zone). The 7- matrices are 
defined as 71 = iay, and 72 = —icx, and ^ „ ^ — ^cr^, 
where cTr^ y ^ are the Pauli matrices acting in sublattice 
space. 

Local interactions in this system can induce collective 
instabilities. Spontaneous symmetry breaking leads to a 
gap in the electronic single-particle excitation spectrum. 
We concentrate on order parameters which break the sub- 
lattice inversion symmetry, i.e., the two sites per unit 
cell become chemically distinguishable. Two possible in- 
stabilities correspond to charge-density wave (CDW) and 
spin-density wave (SDW) states, whose order parameters 
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in terms of electronic operators read 

X = (Xs,Xt) = ((*<T,K*a.K), a^nfaa'-^ a' , (2) 

where t acts in spin space and is simply a vector com- 
posed of Pauli matrices. In the above representation 
the singlet Xs corresponds to staggered charge density 
(CDW), and the triplet xt to staggered local magneti- 
zation (SDW). Microscopically, a finite Xs ixt) may be 
induced by a large nearest-neighbor (on-site) repulsion. 
In the following we explicitly study the case of CDW 
and comment on differences with the SDW case as we go 
along. 



since it is needed for the generic transport equations. 
It is important to note here that L lacks the Lorentz 
symmetry due to the bosonic velocity Vg ^ vp. 

Interestingly, the above Gross-Neveu theory is renor- 
malizable in 3 -I- 1 dimensions, where both the Yukawa 
coupling gs and the self-interaction coupling Ug become 
dimensionless. It is then possible to study the effective 
Yukawa theory in d = 3 — e dimensions using renormal- 
ization group (RG) methods with e being the control pa- 
rameter for the perturbative series (e = 1 corresponds 
to the two-dimensional model). We note that other au- 
thors have considered similar field-theories using an RG 
framework in the context of d-wave superconductors'^4'. 



A. Model I: Yukawa-theory and the description of 
the semimetal-to-insulator transition in the 
framework of an e-expansion 

We start with introducing model I, which is a field 
theory of Dirac fermions coupled to a Landau- Ginzburg 
order parameter theory via a Yukawa-type coupling. For 
the following discussions we generalize the model to d 
dimensions and adopt a Lagrangian formulation. The 
generalization to d dimensions eventually allows to treat 
the bosonic and fermionic sectors on an equal footing in a 
perturbative renormalization group treatment. We gen- 
eralize the Dirac Hamiltonian to three dimensions such 
that matrix elements have a structure like in two dimen- 
sions: on a technical level this implies that we evalu- 
ate angular integrals as if they were in two dimensions, 
but keep the integrals over absolute values in d dimen- 
sions. The following discussion closely mimics the main 
steps of Ref. [TOl We decompose the local interaction 
into separate dynamical fields via Hubbard-Stratonovich 
decoupling and formulate the effective action near the 
semimetal-to-insulator transition as S" = J drdxL, with 
L = L f+Ly+Lb+Lc-, where the bosonic part is described 
by 
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^''V^+i) Xs+Us{XsXsf 



(3) 



and the Yukawa term that couples bosonic and fermionic 
fields reads 



9sXs 



(4) 



We note that Lf, describes the CDW fluctuation mode, 
where t is the tuning parameter, Us the self-interaction, 
and Vs the velocity. It is thus a Landau-Ginzburg theory 
for an Ising order parameter. The physical meaning of 
the Yukawa term is that a CDW can decay into a pair of 
fermions. We also introduce the inverse bosonic propa- 
gator at this point 



y, r, r') = {-d^ - v^^ + t) 5{x - y)5{T - r') 



(5) 



1. RG description of tfie semimetal-to-insulator transition 

In the following discussion we sketch the RG treatment 
of Ref. JO, We first restore Lorentz invariance of the 
model, implying we set — vp = 1 throughout this 
paragraph. This is a crucial point and we extensively 
comment on this step below. 

We take L as the starting point and integrate out 
fermionic and bosonic modes within the momentum- 
shell defined by A/6 < (w^ ^2^/2 ^ ^_ further- 



more introduce dimensionless couplings g = 
with e = 3 — d yielding 



— 9s 



and 



dl_ 

dhih 

du 
d\nb 



g^ie-lg'^) and 
u(e-8.g2)_ 36^2 + 2/ , 



(6) 



to one-loop order. The above flow equations obtain for 
the bosonic mass t tuned to criticality. The intricate 
flow diagram is schematically depicted in Fig. [2] Most 
interestingly, the usual Wilson-Fisher flxed point (WF) of 
the (/)^-theory is unstable towards finite Yukawa coupling 
g leading to an additional non-trivial fixed point called 
SMI located at 



and 



63 



(7) 



This fixed point describes the semimetal-to-insulator 
transition and the associated tuning parameter is the 
bosonic mass parameter t (Eq. ([s])). The correct the- 
ory thus is described by electrons and bosons. A simple 
order parameter theory would not suffice to describe the 
physics correctly, since there the critical point would es- 
sentially be given by WF. 

In the above picture we have assumed the equivalence 
of bosonic and fermionic velocities, i.e., Lorentz invari- 
ance. A closer look at the theory reveals that this is 
only asymptotically correct and breaking the relativistic 
invariance actually corresponds to an irrelevant pertur- 
bation with respect to the fixed point SMI. The corre- 
sponding RG equation was derived in Ref. 1Q_ and reads 



dd 
dlnb 



4e 



(8) 
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FIG. 2: Schematic RG flow diagram: g and u are the Yukawa 
and the bosonic couphng, respectively; G denotes the un- 
stable Gaussian fixed point of the Landau theory and WF 
the Wilson-Fisher fixed point, which is unstable towards SMI 
(semimetal-to-insulator fixed point) upon switching on the 
Yukawa coupling g. The transition is tuned by the bosonic 
mass parameter t which is perpendicular to the u — g plane 
and not shown explicitly. 



at 5*, where 5 — 



This parameter will play a vital 



role in the analysis of the leading scattering mechanism 
in Sec. IIIIA2I 

We note that in the case of SDW the RG equations 
remain structurally intact, with modified combinatorial 
factors coming from the order parameter being 0(3) in- 
stead of Z2. 



2. Long-range tail of Coulomb interaction 

It is know that in a half-filled electronic system on 
the honeycomb lattice long-range Coulomb interaction 
is marginally irrelevanlP^. It is thus also important to 
understand the role of long-range Coulomb interaction 
at SMI. In order to do this one has to generalize the 
Coulomb interaction to d dimension following Ref. I9I10I 
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(9) 



The physics of the quantum-critical point SMI has 
been analyzed in great detail by Herbut et and we 
only repeat one key feature. A prediction of the field the- 
ory is that right at the quantum-critical point the fermion 
propagator behaves as 



GJ^ cx {u? + k^y 



(11) 



where r]f = yf 11^6 anomalous dimension. In this 
sense there are no well-defined fcrmionic quasiparticles 
at the critical point. However, it was shown by Damle 
and Sachdev^"'^ that a Boltzmann transport analysis in 
terms of quasiparticles still captures the relevant physics 
in the hydrodynamic regime {huj <§; ksT). In that sense 
the semiclassical Boltzmann equation only constitutes 
a convenient starting point to solve the more compli- 
cated Quantum-Boltzmann equation, which we refrain 
from showing in full bloom in this paper. 



B. Model II: the Gross-Neveu model and 
description of the semimetal-to-insulator transition 
in large-N 

In the following we consider an alternative formulation 
of the physics described above. The advantage of the be- 
low model is that we can analyze the model directly in 
two dimensions. In order to do so we employ a general- 
ization in the spirit of the large-N method. The large-N 
limit is obtained by generalizing the fermionic spinors 
with a new index i = l,...,iV. The model of the free 
electrons was introduced before in Eq. ([ij . In the formu- 
lation of the model we follow Ref. [H The interaction is 
now modeled by contact interaction of the type 

Lint = ^ a,K,i) 



(12) 



The physical equivalence of the above model with model 
I is established upon noting that gs drives an instability 
towards CDW, whereas gt drives the instability towards 
SDW. 



This form ensures that integration over ao induces a oc ^ 
density-density interaction between fermions in any di- 
mension. It is found^" that Coulomb interaction obeys 
the following RG equation 



dlnb 



327r^ 



l)e^ 



(10) 



One observes that our choice of Lc (Eq. ^) indeed pro- 
duces a marginally irrelevant long-range interaction in 
any dimension. This establishes long-range Coulomb in- 
teraction as the least irrelevant coupling at SMI. 



1. Quantum- criticality on the level of large-N 

In the following we again concentrate on the case of 
the CDW, implying that we set gt — 0. For large N, gs 
can give rise to a dynamically generated mass of the type 



which is determined by the gap equation 



9s 



duj 



^'-"^ d^k 



(13) 



(14) 
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which can be integrated to yield 



1 = 




(15) 



where A 3> to is the ultraviolet cutoff. Demanding in- 



variance of m under rescaling A 4 we obtain 



dg 
d\nb 



= -g-9 



(16) 



where 



2g.A 



g was used (note that Eq. ( 15 ) differs 



from Ref. IHldue to a slightly different choice in the cutoff 
scheme), g* = —1 represents a quantum-critical point, 
where for g < —1 runaway flow into a phase with broken 
chiral symmetry, i.e., m 7^ obtains. This model repre- 
sents an alternative starting point for studying the crit- 
ical transport properties in the framework of the Boltz- 
mann equation. 



2. The role of long-range Coulomb interaction 

Without going into the details of the derivatiorP we 
simply state that the effect of long-range Coulomb inter- 
action has also been studied in the context of this model 
in the limit of large-N. In order to do so the Coulomb 
interaction (|9| for c? = 2 is equally generalized to large 
N and the flow equation of the charge reads 



d\nb 



2nN 



(17) 



which confirms that the long-range Coulomb interaction 
is marginally irrelevant, in qualitative agreement with 
Eq. jTol) from model I and Ref. IH 



We close this section with a generic comment about re- 
cent numerical simulations performed for Dirac fermions 
with only long-range interactions in terms of a lattice 
gauge theory simulation. It was showrP^ that even with- 
out short-range interactions a quantum phase transitions 
as a function of the fine-structure constant can be driven 
from a semimetal to an insulator, where in the insulat- 
ing phase chiral symmetry is observed to be broken. This 
transition happens at a critical value of the fine structure 
constant and numerical values for this critical coupling 
come close to the value expected in vacuum. Within 
numerical accuracy the transition proves to be second 
order and more remarkably the critical exponents seem 
compatible with the above two models with short-range 
interactions exclusively. 

We note that for a purely Hubbard type model there 
has been a recent work discussing the Mott transition 
which is not driven by collective symmetry breaking, but 
a transition into an intermediate paramagnetic spin liq- 
uid phase is observecP^. Only upon increasing the local 
interaction the model has a transition into an antifer- 
romagnetic phase. The nature of this spin liquid phase 



is currently not fully understood and not subject of our 
discussion. 



III. BOLTZMANN TRANSPORT EQUATION 

The Boltzmann transport equation in combination 
with RG methods and large-N has proven a valuable tool 
in a nalyzing transport properties of quantum-critical sys- 
tem422tl23|25j j^^ |.]^jg section we introduce the basic nota- 
tion and concepts for the calculation of the critical con- 
ductivity following closely Ref. [55l 

It is well known that in momentum space one can di- 
agonalize the Hamiltonian ([T]) using a unitary matrix of 
the form 



K* -K* 



with 



(18) 



(19) 



and k = |k| = \K\. Expressing the Hamiltonian Hq in 
terms of particle and hole operators 



7+,<T,«:(k) 

7-.<T,K(k) 



f>k'*..K(k) 



(20) 



we obtain the diagonal form 



d^k 

(2^ 



rAfc 



7l,.,Jk)7,,,,,(k), (21) 



where the sum A extends over electron- and hole-band, 
denoted — . 

We can also express the electrical current in terms of 
the 7± . For the case of a spatially independent current 
the result can be written as 



J = J/ + 3 II 



with 



d^k Ak 

(27r)2T 



(22) 



(23) 



J/ measures the current carried by the motion of the 
quasiparticles and quasiholes; the prefactor A accounts 
for these excitations having opposite charges. The oper- 
ator 3 II is not shown explicitly; it describes the creation 
of a quasiparticle-quasihole pair (it corresponds to the 
so-called Zitterbewegung, see Ref. [T]) 

As in the problems studied in Refs. [20U23|25t in a 
particle-hole symmetric situation a current carrying state 
with holes and electrons moving in opposite directions 
has a vanishing total momentum and the current can 
decay by creation or annihilation of particle-hole pairs 
without violating momentum conservation. This is the 
physical reason why at the particle-hole symmetric point. 
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i.e., at vanishing deviation of the chemical potential from 
the Dirac point, the d.c. conductivity is finite even in 
the absence of momentum relaxing impurities. However, 
at finite deviation from particle-hole symmetry, driven 
by chemical potential or the presence of second neigh- 
bor hopping t' , a driving electric field always excites the 
system into a state with finite momentum which cannot 
decay. This entails an infinite d.c. conductivity in ab- 
sence of Umklapp scattering or disorder. 

As a first step one defines the distribution function 

/A(k,t) = (7L,Jk,i)7,^^^k,i)). (24) 

Note, there is no sum over tr, k on the r.h.s.. We fur- 
thermore assume that the distribution functions are the 
same for all valleys and spins, which has to be expected 
for symmetry reasons. In equilibrium, i.e., in the absence 
of external perturbations, the distribution functions are 
Fermi functions 



/+(k,t) = f{k) = 
f.{k,t) = f{-k) 



e(fe-A')/T + 1 ' 
I 

~ g(-fe-A<)/T _^ I 



(25) 



where we temporarily allow for a finite chemical potential 
/i. 

In principle, off-diagonal elements such as (7^-7^) are 
also created by an electric field. However, they are not 
needed to evaluate J/, which is the part of the current 
that we focus on in the hydrodynamic regime, lo <C fc^T. 
Furthermore, the off-diagonal elements feed back to the 
kinetic equation of the diagonal elements only to higher 
order in perturbation theory. In the presence of an exter- 
nal electric field E acting as a driving force on quasielec- 
trons and quasiholes, we find the semiclassical Boltzmann 
equation 



d_ 
Ft 



eE- 



d_ 



/A(k,<)=/< 



coll 



(26) 



where /coii is the scattering integral. We parameterize 
the change in f\ from its equilibrium value bj^^fil 



/A(k,c.) 



= 2TTS{uj)f{Xk) 

k • E{oj) 

+ XevF 



/"(Afc)(I 



f{Xk))gik,Lu) , 
(27) 



where we have performed a Fourier transform in time 
to frequencies, w, and introduced the unknown function 
g(k,ijj). At the particle-hole symmetric point (fi — 0), an 
applied electric field generates deviations in the distribu- 
tion functions having opposite sign for quasiparticles and 
quasiholes. Formally, this is a consequence of the driving 
term in Eq. (26) being asymmetric under A — > —A, and 



thus the solution has to be asymmetric as well. This re- 
flects the fact that there is an increased number of quasi- 
holes and quasiparticles moving parallel and antiparallel 
to the exciting field, respectively. As quasiparticles and 
-holes have opposite charges, their electrical currents are 
equal and add up, while their net momenta have oppos- 
ing signs and subtract to yield zero. In this paper we 
consider instabilities, which do not break particle-hole 
symmetry, and thus the above parametrization is strictly 
justified. 



A. Boltzmann equations for model I 



We can formulate the Boltzmann equation for the elec- 
trons and holes (we only show the equation for the elec- 
trons; for holes it follows by symmetry) to lowest order 
in the electron-boson coupling as 



e^--Edkfik) = Agl j ^Af++(k, k - q) Im I?(q, - |k - q|) x 

X {n,{k - |k - q|) (/+(k) - /+(k - q)) + /+(k) (1 - /+(k - q))) 
+ 45.^1 J^M+_(k,k-q)Imi?(q,fc+|k-q|)x 

X (nb(fc + |k-q|)(/+(k)-/_(k-q)) + /+(k)(I-/_(k-q))) , (28) 

I 



where with 



T,,,(k,k') = [u^'t^U^, 



M,,,(k,k') = r,v(k,k')T,-;i(k,k') (29) " 2 I' kk' ' ^^"^ 
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D is the bosonic propagator ^ and ni,{e) the Bose dis- 
tribution function. 

There exists an analogous equation for the bosonic 
degrees of freedom. Since the bosons do not couple to 



the electromagnetic field, the driving term (l.h.s.) is ex- 
actly zero. However, the scattering integral couples the 
non-equilibrium fermions to the bosons constraining the 
bosonic distribution function to solve 



= 



(27r)° 



MH.+ (k,k + q) K(k)(/+(k + q) + /_(k + q) - (/+(q) + /_(q)))) 



d'^q 

W) 

d'^q 



I (1^^^+-^'''^ + (^ib(k)(/+(k + q) + /-(k + q) - (/+(q) + /_(q)))) 
, Af++ (k, k + q) (/+ (k + q)(l - /+ (q)) + /_ (k + q) (1 - /_ (q))) 



/ (^^'^+-^''' ^ + ^^^^^ + "1^^^ " ^-^"^^^ + /-(k + q)(l - /+(q))) 



r 



(31) 



This bosonic Boltzmann equation (31) has to be solved 
simultaneously with the fermionic one (28 1. In principle. 



there is also a boson only scattering term coming from 
the self-interaction (xf-termjf^S. This term is not shown 
explicitly for reasons of conciseness. It turns out that its 
neglect proves to be exact in the case at hand. This is 
rooted in the fact that the bosons remain in equilibrium, 
which is shown explicitly below. 



1. The effect of 'boson-drag' 

Severe complications in calculating electromagnetic re- 
sponse functions in systems in which fermions and col- 
lective (slave-)bosons are coupled to each other result 
from the fact that one has to solve coupled Boltzmann 
equations (Eq. 



p8| ) and Eq. ([31])). The physical reason 
for that is simple: electrons transfer momentum to the 
bosons, which themselves are driven out of equilibrium. 
Since the bosons also scatter from the fermions, there 
is a feedback effect which in principle has to be taken 
seriously. Often this is circumvented by assuming the 
bosonic sector to be in equilibrium, which is equivalent 
to saying that the bosonic sector equilibrates due to lat- 
tice effects or impurities on a time scale faster than the 
electronic scattering time. This can be an oversimplifying 
assumption with no a priori justification. One can ap- 
preciate this most easily by considering a Galilean invari- 
ant electronic system interacting with effective bosonic 
modes. In this situation neglecting the generic non- 
equilibrium situation for the bosonic sector yields fun- 
damentally wrong results, namely a finite d.c. conduc- 
tivity. One can convince oneself that taking into account 
the Boltzmann equation for the bosons reestablishes the 
physically correct result of infinite d.c. conductivity (see 
chapter on 'phonon-drag' in Ref. rSB) . 



We now turn our attention to Eq. (31) and expand 



order in the applied external field E using Eq. (27). For 
the bosonic sector we stick to the symbolic 



nb(k) =nO(e(k))-f Mk$(fc) 



(32) 



with lik being linear in the applied field and n° (e) as the 
equilibrium Bose distribution function. A simple analysis 
shows that $(fc) drops out of Eq. (31 ) and the remaining 
parts annihilate exactly, as they should for consistency. 
The above scattering integral is thus solved exactly by 
letting 



n,(k) = nO(e(k)) 



(33) 



and using Eq. (27) for the fermions. We have thus shown 



the fermionic and bosonic distribution functions to linear 



that at the Dirac point the 'boson-drag' vanishes exactly 
and considering the bosons in thermodynamic equilib- 
rium is exact. 

One can convince oneself that in a situation with finite 
chemical potential /i this ceases to be true and the bosons 
have to be treated as being out of equilibrium. In such 
a situation a finite amount of disorder is needed to relax 
the net momentum and the associated current. 

We have thus established that in clean graphene at 
fi = the bosonic sector remains in equilibrium. We can 
understand this more intuitively as follows: we first con- 
sider an electronic band with quadratic dispersion inter- 
acting with phonons. In this situation the bosonic sector 
absorbs some momentum from the electrons which scat- 
ter and is itself driven out of equilibrium. In graphene 
we have no quadratic disperion and by virtue of the per- 
fect particle-hole symmetry (/i = 0) the charge carrying 
modes at the Dirac point are modes with net momentum 
zero. Thus the bosons do not absorb any momentum 
which they have to dissipate. This implies they are not 
driven out of equilibrium and treating the bosons in equi- 
librium is no approximation but exact to linear order in 
the applied external field E. Note that the above reason- 
ing is exclusive to electrical transport and does not hold 
for thermal transport. 
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This hne of argument is not exclusive to graphene or 
the matrix elements of the above interaction, but also 
applies to other systems with perfect particle-hole sym- 
metry which is unbroken by the interaction mediated by 
the bosons. As^a timely example one could consider bi- 
layer graphen^ where the quasiparticles interact with 
effective bosonic degrees of freedom and the reasoning 
still holds. 



2. Cancellation to infinite order 



teraction constitutes another perturbation to SMI, which 
is only marginally irrelevant and consequently vanishes 
more slowly as temperature is lowered, namely in a log- 
arithmic manner. From the above analysis we conclude 
that the dominant current relaxation mechanism stems 
from long-range Coulomb interaction. This implies that 
to leading order the behavior in the semimetal is the same 
as at SMI. 

This statement in itself is very surprising. We are now 
going to reexamine and reinterpret it in the framework 
of model II. 



For dimensional reasons one expects the solution to 
Eq. ( 28 1 to yield a quantum-critical conductivity which 
is independent of temperaturd^illSS] turns out that this 
is not true. To quadratic order in the interaction g the 
scattering integral ( 28 1 is exactly zero due to a conspiracy 



of the matrix element and energy conservation. This in 
itself does not rule out a temperature independent con- 
ductivity, since contributions to all orders in g provide 
scattering times independent of temperature. However, 
one can convince oneself that the kinematical constraint 
shown in Fig. |3] annihilates contributions to current re- 
laxation due to the Yukawa coupling g to all orders in 
perturbation theory, provided vp = Vg- This can easily 



be seen by looking at the matrix elements Eq. ( 30 ) and 
imposing energy and momentum conservation as shown 
in Fig. |3] The exact cancellation is only operational for 
Va = Vp. Breaking of Lorentz invariance, i.e., allowing 
for Va ^ Vp allows to relax the current. This implies 
that current relaxation due to the collective CDW de- 
gree of freedom is determined not only by gs but also by 
an irrelevant parameter, namely (5, which was introduced 
as the dimensionless parameter measuring the breaking 
of Lorentz invariance. Since ga assumes its fixed point 
value, which is a pure number, the inverse scattering 
time vanishes upon lowering temperature at least like 
E oc T*/'', see Eq. Q, leading to a diverging d.c. con- 
ductivity. One now has to determine the least irrelevant 
scattering mechanism among all possible ones. As we 
argued before, all scattering events coming from the in- 
teraction with the collective mode scale with additional 
powers of 8 and thus are irrelevant. 




FIG. 3: Momentum and energy structure of the Yukawa ver- 
tex g for scattering of quasiholes (A = — , full line) and quasi- 
electrons (A = — ) from the collective mode (wiggly line). One 
can work out that for Vs — vp the bosonic spectral function is 
composed of sharp quasiparticles and only coUinear scatter- 
ing is allowed; coUinear scattering, however, annihilates the 
associated matrix elements. 



B. Boltzmann equation for model II 

The following analysis is carried out directly in 2 di 
mensions starting from the interaction term Eq. ( 12 1 
The collision term on the r.h.s of Eq. (26) can 



We have seen in Sec. II A 2 that long-range Coulomb in- 



be de- 
termined using Fermi's golden rule or equivalently from 
a Quantum-Boltzmann equation. We refrain from an ex- 
plicit derivation and refer the reader to Ref.'^^, which 
gives an explicit derivation for the case of long-range in- 
teractions which is easily generalized to the present case. 
The generalization leads to 
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where 



dt 



eE.|)/.(k,i) = - 



(27r) r (fki (Pq 
^ J (27r)2 (27r)2 



5{k - fci - |k + q| + |ki - q|)i?i(k, ki, q){/;,(k, t)/_A(ki, - fx{^ + q, i)] 

X [1 - /_A(ki - q,0] - [1 - /A(k,i)][l - /_A(ki,t)]/A(k + q,i)/_A(ki - q,t)} 
^(fc + fci - |k + q| - |ki - q|)i?2(k, ki, q){/A(k, t)/A(ki, i)[l - /A(k + q, t)] 

X [1 - A (ki - q, t)] - [1 - A (k, t)] [1 - A (ki , t)]h (k + q, t) h (ki - q, t) } I , 



i?i(k,ki,q) = 4Ar(|r+__ + (k,ki,q)|' + |r+-+-(k,ki,-k-q + ki) 
i?2(k,ki,q) = 47V|r++++(k,ki,q)|' . 



(34) 



(35) 



Note that the colhsion kernels Ri and i?2 both come with 
a prefactor N. There are additional contributions com- 
ing from crossed diagrams, which are however down by 
a factor 1/A^, see Ref. i25t The terms proportional to Ri 
represent collisions between oppositely charged particles, 
while those proportional to i?2 are collisions between like 
charges. Other processes, where a particle-hole pair is 
created, turn out to have vanishing phase space: it is not 
possible to fulfill momentum and energy conservation at 
the same time. This is a peculiarity of the linear disper- 



sion, i.e., Ek = vpk (see also Ref. ES]) . The scattering 
matrix elements are given by 



TaAi A.A3 (k, ki, q) = |^Taa3 (k + q, k)TA, A. (ki - q, ki) 



(36) 



We now proceed to the linearization of ( 35 1 by inserting 



the parametrization (27) (we define g\{k, uj) = Xg{k,uj)) 
and find 



{-iLogx{k,uj) ~ Xvf/T) k (27r) f(fki (fq 



(e"Ffe/T + i)(e-^'Ffe/T + i) k " 
5{k-ki - |k + q| 



VF J (27r)2(27r)2 
|ki -q|)i?i(k,ki,q) 



ki 



(k + q) ,,, I , (ki 



fci"""^ ' |k + q 
S{k + fci - |k + q| - |ki - q|)fl2(k, ki, q) 



ki 

fci 



, (k-l-q) ,., . , (ki 



ki 



-5_A(|ki - q\,uj) 



-gx{\ki - q\,uj) 



(37) 



For the subsequent discussion we take the limit uj — 
and focus on the solution of the linearized transport equa- 
tion in Eq. ( [37| for the function g. The following discus- 
sion closely follows the presentation given in Ref. 25, We 
can view the right hand side of Eq. ( 37 1 as a linear op- 



erator, the so-called collision operator C, acting on the 
function {k/k)g{k). A key property of C is that it is 



Hermitian with respect to the natural inner product 
(51I52) = Xl/ (^jj^5i,a(%2,a(^) ■ 



(38) 



This Hermiticity followJ^^ from symmetry properties of 
Ri and R2 under exchanges between incoming and out- 
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going momenta, which are very similar to those used in with the exphcit form of the functional given by 
establishing Boltzmann's H-theorem. 

Related to the above properties of the collision op- 
erator, we can introduce a functional Q[g], such that 



Eq. (37) is equivalent to finding its stationary point 



(39) 



8vF J (27r)2 (27r)2 (27r)2 | 

6ik - fci - |k + q| + |ki - q|)fii(k, ki, q) 

(g-i>Ffc/T _^ i)(e''-Ffci/r + l)(e^F|k+q|/r _|_ i^dg-i-Flki-ql/T _^ 

k ,,^ ki ^ (k + q) ,1, I, (ki - q) 

pik) - ^5(^1) - IkT^^dk + ql) + - q 

Sjk + fci - |k + q| - |ki - q|) j?2(k, ki, q) 

k ki ^ (k + q) ,1, I, (ki - q) 

p{k) + ^5(^1) - IkT^-^dk + ql) - - q 



d2fc 



g{k)vF/T 



(27r)2 (e«Ffc/T 4, i)(g-^pfc/T + 1) ■ 
I 



(40) 



In the previous analysis of a quantum Boltzmann equa- 
tion for massless Dirac fermions in two dimension&22 , it 
was noted that the phase space for scattering of particles 
was logarithmically divergent in the collinear limit. For 
the interaction considered in that paper, the collinear 
scattering cross-section vanished, and so this singular 
phase space density had no important consequences. The 
collinear scattering does not vanish for the present inter- 
action, and so we need to take this logarithmic divergence 
seriousljl^SlMI. 

The physical origin of the divergent collinear scat- 
tering is related to the linear dispersion which implies 
that quasiparticles or -holes moving in the same direc- 
tion share the same group velocity, independent of their 
energies. This leads to a diverging duration of colli- 
sions of nearly collinear particles, which is enhanced due 
to the low space dimensionality. To the extent that 
collinear scattering is very strong, and considering fre- 
quencies much smaller than the inelastic scattering rate, 
we may expect that quasiparticles and holes that move 
in the same direction in the plane will establish a pseudo- 
equilibrium characterized by an effective chemical po- 
tential and an effective temperature which, however, de- 
pends on the direction of motion. 

In linear response the deviations of these effective pa- 
rameters from the equilibrium values and T have to 
vary with k/fc • E for symmetry reasons. The remaining 
dominant mode of the function g will correspond to an 



effective shift in chemical potential which translates into 

(41) 



where the prefactor has been chosen so as to make xi^) 
dimensionless. With this Ansatz, which will be confirmed 
below, it simply remains to determine the prefactor x, 
yielding the leading term in the non- equilibrium distri- 
bution. The whole reasoning presented here is given in 
much more detail in Refs. [25>28i Note that the effective 
chemical potential shift ranges between ±xfiVFeE/T de- 
pending on the direction of motion. Comparing this to 
the temperature allows us to estimate the threshold elec- 
tric field strength, eEwa — T^/hvp, below which non- 
linear effects should remain small. 



The Ansatz (41) can be justified along the lines of 
Ref. [251 The phase space for scattering diverges upon 
the perpendicular momentum going to zero. However, 
it turns out that the scattering matrix is not zero for 
collinear scattering. This implies that the logarithmic di- 
vergence for the perpendicular momentum going to zero 
has to be cut off by higher-order self-energy corrections 
to the free fermions. This implies that modes which are 
not zero for the collinear scattering channel have an ad- 
ditional factor In j- , which in the limit ^ 1 is a large 
number. 

It thus turns out that the diverging phase space for 
collinear scattering helps to essentially solve the problem 
exactly in the limit ln(A^/gs) ^ 1. 
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We conclude that up to corrections of order 
[\n{N/ gs)]~^ , we can choose g to be of the form 



9ik) 



rp2 



C. 



(42) 



We insert this parameterization into the functional Q[g] 
in Eq. (401; the solution of the stationarity condition in 



Eq. ( 39 ) is then equivalent to requiring the vanishing of 



the derivative with respect to C. We numerically evalu- 
ated the integrals in Eq. ( 40 ) using an elliptic co-ordinate 
system to solve the energy conservation constraint!^, and 
obtained 



S[.9] 



1 In 2 

T 4^ 



- 2C 



with K^l. 8448 x 10"'' ± 1.0 x 10" 
arity condition we then obtain 



C 



nglT'^ 



From the station- 



(43) 



The conductivity can be obtained from C by combining 
Eqs. ([23l), (pTl and 



ct(w 0) - 



(44) 



where we have re- inserted factors of h and ks- This 
seems to diverge upon lowering temperature, however, 
we have to keep in mind that at criticality we have 
g* = = —1- Identifying the running energy scale 

A with the temperature fc^T we end up with 



a{Lo = 0) 



AN^ In 2 

h KTT^ 



(45) 



This result is independent of temperature, as we expect 
for a quantum-critical point in two dimensions. It is 
worthwhile noting that the prefactor is oc N'^, which is 
different from the large-N transport analysis of the 0{N) 
rotor model shown in Ref. [201 which only has a prefactor 
of N. The difference stems from the fact that in our case 
all N flavors couple to the electromagnetic field, whereas 
in the latter case only two components were coupled to 
the electromagnetic field. So in contrast to the analy- 
sis of the transport equations in model I we find that 
the universal conductivity is finite, however with a nu- 
merically huge prefactor. We will now comment on the 
relationship between the two results. 



C. Comparison of results in model I and model II 



point, in model II we find a finite scattering due to the 
electron-electron interaction. We will argue in the follow- 
ing that this can actually be interpreted consistently. In 
order to do so we have to analyze the reason for the kine- 
matical constraint in model I. This discussion can most 
easily be carried out directly in two dimensions, where 
for the moment we put all subtleties about controlling 
the perturbation series aside. If we consider the r.h.s 



of Eq. ( 28 1 we realize that the reason for knocking out 



scattering processes stems from the fact that the order 
parameter fluctuations are described by sharp quasiparti- 
cles, i.e., the spectral function assumes the form of delta 
peaks. One can now ask in which sense this describes 
the right physics. In order to do so we can pretend to 
derive model I from model II via Hubbard-Stratonovich 
transformation and integrating out electrons. It turns 
out that the propagator of the order parameter field is 
actually not given by 



L'(z^„,k) oc 



1 



k2 



but instead by 



£)(^'„,k) oc 



(46) 



(47) 



This has severe consequences: while Eq. (46) implies 



that the bosonic spectral function is sharply peaked. 



Eq. ( 47 1 describes a continuum, which still has its maxi- 



mum weight at the former resonance. While a lot of the 
weight still is kinematically blocked, using Eq. (47 1 in the 



Boltzmann approach would immediately render the r.h.s. 
of Eq. (28) finite, albeit small. In that sense the bosonic 



sector in model I does not faithfully describe the charac- 
ter of the bosonic spectral function in the real two dimen- 
sional system. This implies that in reality we expect the 
real bosonic spectral function to still be peaked around 
the resonances, but also to have a continuum background. 
The Boltzmann analysis in model II backs up this state- 
ment. Directly in two dimensions we find a finite temper- 
ature independent inverse scattering time from interac- 
tions with the collective field, albeit with a tiny numerical 
prefactor. The prefactor being so small can consistently 
be interpreted as a remnant of the complete kinemati- 
cal blocking in model I. Our conclusion thus is that the 
Boltzmann equation as it is presented for model I simply 
overestimates the kinematical blocking by annihilating 
the scattering completely. We thus conclude this section 
by stating that the asymptotic quantum-critical univer- 
sal conductivity is numerically huge, but independent of 
temperature, as we expected to find. 



We have seen that model I and model II seemingly 
give contradictory answers to the question whether there 
is a finite universal conductivity. Whereas in model I we 
find that a kinematical constraint knocks out scattering 
from the order parameter fiuctuations right at the critical 



IV. QUANTUM-CRITICAL CONDUCTIVITY 

We will now discuss the implications of the results of 
Sec. |III[ We choose to carry out the discussion in the 
framework of model II for reasons of enhanced clarity. 
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We found that the universal quantum-critical conductiv- 
ity within model II is a pure number times /h. This 
central result is shown in Eq. (45 1. However, as we ar- 



gued before, the prefactor of the inverse scattering time 
is numerically tiny, consistent with the complete blocking 
found in model I. This naturally leads us to consider the 
influence of the remaining marginally irrelevant interac- 
tion, which is given by the long-range Coulomb interac- 
tion. 

We established that in the clean system at zero chemi- 
cal potential there are two main sources of current relax- 
ation namely (i) short-range interaction which drives the 
phase transition and (ii) long-range Coulomb interaction. 

In the following we compare the inverse scattering 
times due to the two independent processes. The dis- 
cussion is simplified by the fact that we can restrict our 
discussion of the inverse scattering times to one hydrody- 
namic mode, see Sec. |III[ We know that in the framework 
of the Boltzmann equation collision integrals add. Thus, 
by virtue of the most relevant relaxation processes for 
both types of interaction living on the same hydrody- 
namic mode inverse scattering times also simply add up. 
This is what is usually called Matthiessen's rul#^ within 
the relaxation time approximation, but here it is essen- 
tially exact. In the following we again lend from Ref. l25l 
and simply take the numerical value of the inverse scat- 
tering time on the zero mode and compare it to the one 
obtained in Eq. ( 45 1 . Since for practical purposes the 



large-N parameter is a number of order on we can 
safely drop it from our discussion in the physically rele- 
vant situation. Schematically, we can write the collision 
integral composed of the sum of the two processes as 



1 



= {a + ba^{T)) T 



where 



a{T) = 



1 



^In^ 

4 A 



(48) 



(49) 



is the running long-range coupling. The perturbation 
theory of both terms is formally controlled in the small 
parameter 1/A^. However, there is a big difference: while 
the first term in brackets is temperature-independent, the 
second term scales to zero logarithmically upon lower- 
ing temperature. This implies that asymptotically in the 
very low temperature limit the second term vanishes and 
the first term determines the universal temperature in- 
dependent quantum-critical conductivity, in agreement 
with the previous discussion. For higher temperatures, 
however, the contribution due to long-range Coulomb in- 
teraction will win. 

We can make a very crude order of magnitude esti- 
mate for the temperature T^r, where the crossover from 
universal quantum-critical conductivity to the regime, in 
which conductivity is determined by long-range Coulomb 
interaction, takes place. We know that in graphene the 
ultraviolet cutoff is of the order eV, which translates to 
A w IG^K. Furthermore, the bare long-range Coulomb 



interaction is a di mensionless number of order one. We 
showed in Sec.|m]that a w 10"^ and Ref.|25]found 6 w 1. 
We can estimate the crossover temperature Td- as the 
temperature, where both contributions are comparable 
in size. This crossover temperature can conservatively be 
estimated to be of the order ~ 10^^*^/^. This implies 
that for all practical purposes the universal quantum- 
critical conductivity will be masked by the long-range 
Coulomb interaction. 

Interestingly, this also implies that the quantum- 
critical conductivity at the semimetal-to-insulator criti- 
cal point has the same characteristic transport properties 
as the one in the semimetal, i.e., for a gas of hot Dirac 
fermions interacting via long-range Coulomb interaction. 
Furthermore, this implies that the conclusions of the hy- 
drody namic t ransport relations obtained in a series of 
paper^^mUMl g^jgQ hold at the quantum-critical point. 

We simply cite the result for the Coulomb interac- 
tion limited minimal conductivity which has been an- 
alyzed elsewher cj^^ l ^^tt SO] xhere, a Boltzmann equation 
was solved directly in two dimensions. This procedure 
has been presented at length in Ref.j^Sland we only high- 
light the final result for the d.c. conductivity reading 



cr = 0.76 



1 



h a{TY 



(50) 



Using Eq. ( 49 ) in Eq. ( 50 ) we observe that upon lowering 
the temperature the conductivity diverges as a scales to 
zero logarithmically. In Eq. ( 49 1 ao is a number of order 



one and A constitutes the high energy cutoff set by the 
lattice. 

We end this discussion by saying that we expect this 
to be the correct minimal conductivity in region II in 
Fig. [T] for all physically relevant situations. This implies 
that within our approach for realistic temperatures we 
cannot distinguish the minimal conductivity in region I 
from region II in Fig. [T] 



V. DISCUSSION AND CONCLUSION 

In this paper we have investigated quantum-critical 
transport of electrons on the honeycomb lattice in the 
vicinity of the semimetal-to-insulator transition. On a 
technical level, we have used a combination of the semi- 
classical Boltzmann transport equation with RG and 
large-N to analyze two related Gross-Neveu type field 
theories. One of the central results of this paper is given 
in Eq. ( 50 ) which is the conductivity in the quantum- 
critical region (region II in Fig. [T]) for physically real- 
istic temperatures. The most important characteristic 
of this minimal d.c. conductivity is that it is set by 
scattering off the long-range Coulomb interaction, which 
is a marginally irrelevant coupling at the semimetal- 
to-insulator transition. The result is surprising since 
the critical conductivity weakly but explicitly depends 
upon temperature through the flow of the dimensionless 
Coulomb interaction parameter (49). In general, for a 
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Lorentz invariant quantum-critical point in two dimen- 
sions described by an interacting fixed point one would 
expect a conductivity independent of temperature. A 
further remarkable consequence of our calculation is that 
to leading order in temperature one cannot distinguish 
the conductivity in Fig. [ijin region I from region II. We 
furthermore showed that below a very low temperature 
scale Tcr a temperature independent universal conductiv- 
ity takes over, in agreement with expecations. 

We finally go back to Fig. [T] and analyze possible 
crossover scenarios. We know that the quantum phase 
transition microscopically is driven by the interaction pa- 
rameter U . Within the two models we consider this trans- 
lates into different quantities. Whereas in model I the 
quantum phase transition is driven by the bosonic mass 
parameter t in the case of model II the phase transition 
is driven by the interaction parameter g. The crossovers 
of the system are shown in Fig. [l] by doted lines with the 
associated crossover temperature T*{t), which is set by 
T*it) (X t"^ = 

One conclusion of the preceding discussion is that the 
universal conductivity as a function of temperature in 
regions I and II is given by 



cr(r)iij cx 



1 



h a^{T) 



(51) 



Another more formal interesting result of our analysis 
of model I is that it is possible to show that the so-called 
'boson-drag' in our setup is exactly zero meaning that an 
equilibrium treatment of the bosonic degrees of freedom 
within the fermionic collision integral is justified to lead- 
ing order in the applied external field E. This statement 
is more general than that and can be extended to sys- 
tems with perfect particle-hole symmetry, in which the 
collective bosonic modes do not break the particle-hole 
symmetry. 



A. Crossover: quantum-critical-to-CDW-insulator 

For the case of the CDW we have an Ising type transi- 
tion. This implies that due to the absence of Goldstone 
modes symmetry breaking is possible at finite temper- 
atures and a finite electronic gap is stabilized. For the 
sake of simplicity we concentrate on a crossover function, 
which captures the essential physics deep in the quantum- 
critical regime and in the insulating regime. One such 
function is of the type 



a{T) 



CDW 
II^III 



1 



Anie T 



(52) 



with Am being a number. 

The activated form of the conductivity in the insulat- 
ing regime III accounts for the fact that we have a finite 
gap separating the conduction band from a valence band, 
link in an ordinary semiconductor. A typical crossover 
curve is shown in Fig. [4j where a comparison of the 




FIG. 4; Resistivity for regions II (quantum-critical region 
)and crossover from quantum-critical to CDW called II-III 
(CDW) as a function of temperature in units of the resistiv- 
ity at T = Tq. The curve for the crossover II-III stays mostly 
flat over a wide temperature range and eventually takes ofi' 
exponentially. In the quantum-critical region II the resistiv- 
ity approaches zero upon lowering temperature in logarithmic 
manner. The plotted curves are shown for numerical values 
ao = 1 (Eq. ra), Am = 0.2, and t = 0.07. 



quantum-critical resistivity (inverse conductivity right at 
the critical coupling) is compared with the crossover func- 
tion for the crossover from region II to region III. The full 
line in Fig. |4] shows the quantum-critical conductivity, 
whereas the dashed line shows the crossover function. 



B. Crossover: quantum-critical-to-SDW-insulator 

In the case of a CDW instability a discrete Ising de- 
gree of freedom is spontaneously broken and long-range 
order at finite temperature is possible. For the SDW 
an additional breaking of spin rotation symmetry is re- 
quired, thus there can be no true long-range order at 
finite temperature by virtue of the Hohenberg-Mermin- 
Wagner theorem. This implies that the fermionic quasi- 
particles are only truly gapped at zero temperature. The 
question is whether the conductivity can still show signa- 
tures of a pseudogap in the conductivity which diverges 
as temperature is lowered. The role of a finite temper- 
ature pseudogap in scaling functions has been analyzed 
before in the context of the Hertz-Millis theory in two 
dimensional systems'^^. In our qualitative consideration, 
however, we go along the lines of Lee et al\^, who con- 
sidered fluctuation effects at the Peierls transition in one 
dimensional systems. It is most convenient to think of 
the problem at very low temperatures in terms of a spin- 
fermion model like model I. We concentrate in the fol- 
lowing on the low-temperature region which is usually 
called the renormalized classical regim^^Sl, There, the 
self energy of the electrons is given by the zero-frequency 
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FIG. 5: Resistivity for regions II (quantum-critical region 
)and crossover from quantum-critical to SDW called II-III 
(SDW) as a function of temperature in units of the resistivity 
at T = Tq. The curve for the crossover II-III stays mostly 
flat over a wide temperature range and eventually takes off 
as a power-law. In the quantum-critical region II the resistiv- 
ity approaches zero upon lowering temperature in logarithmic 
manner. The plotted curves are shown for numerical values 
ao = 1 (Eq. @) and Am = 0.1. 



Matsubara mode of the bosonic propagator, i.e., 

g(q,t/„) ^ ^2 I ^-2 (53) 

where ^ is the finite correlation length which diverges 
upon lowering temperature. In the framework of a Boltz- 
mann theory this implies that the former inelastic scat- 
tering of the electrons off the bosonic mode has become 
elastic and corresponds to the scattering off long-range 
impurities of the ^-type where r is the distance. In order 

for this to be true we assume %jr- = here, which is a 
non-singular limit in the Boltzmann-equation and corre- 
sponds to the "renormalized classical" regime. The tem- 
perature behavior of such a type of impurities has been 
analyzed befor^^ and it was found that the associated 
transport scattering time behaved as 

7 « ^ • (54) 

Again, we assume a crossover function which correctly 
captures the transport characteristics in the quantum- 
critical and the renormalized classical regime of the form 

The respective crossover function is plotted and com- 
pared to the quantum-critical resistivity with notation 
as above in Fig. [5] 



C. Role of particle-hole symmetry breaking and 
disorder 

One element we have not yet discussed is the role 
of particle-hole symmetry breaking and related to this 
disorder. If the system does not have perfect particle- 
hole symmetry some rigorous statements, such as the 
one about vanishing boson-drag and finite conductivity 
without disorder, cease to be true. In practice, there 
are several ways to break particle-hole symmetry. One 
is to allow for a finite chemical potential, which however 
brings us away from half-filling. Another way is to take 
into account second neighbor hopping t' . If particle-hole 
symmetry is absent, there is no current relaxation in ab- 
sence of impurities, since the momentum mode is excited 
in that case. Also, in that case the boson-drag cannot be 
neglected any more. However, we still expect our theory 
to be valid in some limit. The regime of validity of our 
theory is identical to the one which has been worked out 
in great detail Ref. [28l For a finite chemical potential we 
need /i <C T for our theory to hold. For particle-hole sym- 
metry breaking due to next-nearest-neighbor hoppings, 
such as provided by we have to require that for the 
temperature window we are considering that we are still 
within energies, where the spectrum behaves linearly. 



D. Experimental relevance 

Within this paper we have analyzed the resistivity in 
the vicinity of the semimetal-to-insulator transition on 
the honeycomb lattice. More specifically, we have studied 
the crossover from region II to region 111 in Fig. [4] and 
Fig. [5] 

We now argue that our findings bear some relevance 
for the understanding of experiments carried out on sus- 
pended graphene samples. Remarkably, the curves shown 
m Fig. |4] and Fig. [5] mimic the almost constant behavior of 
the resistivity upon decreasing temperature over a large 
temperature window followed by an upturn reminiscent 
of an insulator as observed by Bolotin et. aP^ on sus- 
pended samples (see Ref. .32, Fig. 4 inset). The samples 
considered in those experiments exhibit ultrahigh mo- 
bilities and enormous mean free paths with respect to 
impurity scattering (concerning impurity scattering the 
samples are ballistic with an associated mean free path 
longer than the sample length). These high mobilities 
and long elastic mean free paths are obtained after re- 
peated annealing removing impurities from the samples. 
This makes the inelastic scattering mechanism discussed 
in this note a viable candidate to set the limiting scatter- 
ing process for the conductivity. Most strikingly, Bolotin 
et al. find that upon lowering temperature the resistivity 
increases. It is tempting to attribute this behavior to the 
existence of a quasiparticle gap. 

It is well known that increasing the interaction on the 
honeycomb lattice, both short- or long-range, can lead 
to a finite quasiparticle gap in the system. Samples on 
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substrates do not seem to show this tendency, but one 
has to keep in mind that suspending the sample mod- 
ifies the dielectric environment and renders the system 
more strongly interacting. One can thus imagine that 
suspending graphene samples modifies the dielectric en- 
vironment strongly enough to drive the system towards 
strong interactions rendering the charge gap finite. 

At finite densities {fj, ^ 0) a qualitative explanation 
of the experimental findings has been given in Ref. [iOl 
However, here we concentrate on the zero doping case 
and to the best of my knowledge the data is not fully 
understood. 

Recently, the data at charge neutrality have been thor- 
oughly reanalyzed by Drut et alW^ who speculated about 
an underlying, possibly interaction driven gap in the sys- 
tem. The authors determined a couple of phenomeno- 
logical parameters to fit the data. More specifically, a 
phenomenological form of the conductivity 

a ^ag + abg (56) 

where the subscript q stands for quasiparticle and bg for 
background was introduced. The background conductiv- 
ity was set by an unspecified scattering process. Thus 
we do not speculate about the nature of this part but 
concentrate on the quasiparticle part. 

For the interpretation and the fitting of abg the authors 
considered a gap as well as another unspecific scattering 
process: concerning the unspecified scattering process, 
the ratio of the inverse scattering time to temperature, 
i.e.j ^ was found to be independent of temperature. Up 
to logarithmic accuracy this is indeed the case in our sce- 
nario, due to the inverse scattering time due to inelastic 
scattering being given by - oc a^(T)T. Furthermore, the 
authors find that within the relevant temperature range 
the mean free path associated with this scattering pro- 
cess is on the order of a tenth of micrometers. This is 
also in agreement with our results. The mean free path 
corresponding to inelastic scattering due to long-range 
Coulomb interaction has been estimated before in Ref.ESl 
and is given by lee ~ ai^)2T[K] ■ sense long-range 

Coulomb interaction is a natural candidate for one of 



the current relaxation mechanisms needed to explain the 
measurements faithfully. 

Concerning the gap of the quasiparticle, the authors 
argued that the data is not only compatible with a hard 
gap but also with a power-law type contribution. The 
hard gap naturally emerges from the CDW scenario, 
whereas the pseudogap behavior stems from the SDW 
scenario. Both mechanisms result in a diverging resistiv- 
ity at low temperatures. 

Obviously, there are additional other scattering 
sources, which could play a role. It turns out that 
phonons due to the high Fermi velocity play a role at 
rather elevated temperatures (T > 150K) and thus are 
not a candidate to modify the above picture of the trans- 
port characteristicPi^. Another candidate would be the 
corrugation or ripples which are observed in suspended 
graphene. It turns out that on a formal level ripples 
share some characteristics with long-range Coulomb in- 
teraction. To our knowledge no thorough study of the ef- 
fect of ripples on the conductivity in suspended graphene 
exists to date and a detailed analysis is called foi^D. 

We thus conclude this discussion by stating that a 
proximity of free standing graphene to a quantum-critical 
point for a transition from the semimetal to an insulator 
described in this note could provide a natural explanation 
for the (temperature) behavior of the resistivity observed 
in the experiment of Bolotin et alW^ on ultraclean high- 
mobility samples. This, however, would also imply that 
signatures of collective symmetry breaking of the SDW or 
CDW type should be observable in those samples, which 
to my knowledge has not been attempted. This provides 
a necessary independent check of the scenario of graphene 
in vacuum being an excitonic insulator. 
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